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Abstract — We explore several reduced-dimension multiuser 
detection (RD-MUD) structures that significantly decrease the 
number of required correlation branches at the receiver front- 
end, while still achieving performance similar to that of the con- 
ventional matched-filter (MF) bank. RD-MUD exploits the fact 
that the number of active users is typically small relative to the 
total number of users in the system and relies on ideas of analog 
compressed sensing to reduce the number of correlators. We 
first develop a general framework for both linear and nonlinear 
RD-MUD detectors. We then present theoretical performance 
analysis for two specific detectors: the linear reduced-dimension 
decorrelating (RDD) detector, which combines subspace projec- 
tion and thresholding to determine active users and sign detection 
for data recovery, and the nonlinear reduced-dimension decision- 
feedback (RDDF) detector, which combines decision-feedback 
orthogonal matching pursuit for active user detection and sign 
detection for data recovery. The theoretical performance results 
for both detectors are validated via numerical simulations. 

I. Introduction 

Multiuser detection (MUD) [ ] is a classical problem in 
multiuser communications, where a number of users commu- 
nicate simultaneously with a given receiver by modulating 
information symbols onto their unique signature waveforms. 
The received signal consists of a noisy version of the su- 
perposition of the transmitted waveforms, and the receiver 
has to detect the symbols of all users simultaneously. While 
there has been a large body of work developed for the 
multiuser detection problem, it is not yet widely implemented 
in practice, largely due to its complexity and high-precision 
AID requirement. The complexity of MUD arises both in the 
analog circuitry for decorrelation as well the digital signal 
processing for data detection of each user. We characterize 
the decorrelation complexity by the number of correlators used 
and the data detection complexity by the complexity-per-bit 
[1]. 

The conventional MUD detection structure consists of a 
matched-filter (MF) bank front-end followed by a linear or 
nonlinear digital detector. The MF-bank front-end is a set 
of correlators, each correlating the received signal with the 
signature waveform of a different user. Hence the conventional 
MUD requires the number of correlators to be equal to the 
number of users. To recover user data from the MF-bank 
output, various digital detectors have been developed. The 
optimal MUD is the maximum likelihood sequence estimator 
(MLSE) [i], which minimizes the probability of error for 
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symbol detection, but its complexity-per-bit is exponential 
in the number of users when the signature waveforms are 
nonorthogonal. The nonlinear decision feedback (DF) de- 
tector [ ] is a good compromise between complexity and 
performance among all nonlinear and linear MUD techniques 
[ ]. This technique detects users iteratively and subtracts 
the strongest user in each iteration. Both the MLSE and 
the DF detectors are nonlinear methods. Linear detection 
requires lower complexity but with a commensurate reduction 
in performance. This technique applies a linear transform to 
the receiver front-end output and then detects each symbol 
separately. Linear MUD techniques include the single-user 
detector, the decorrelating detector and the minimum mean- 
square-error (MMSE) detector [ ]. When the user signature 
waveforms are correlated, the performance of the single-user 
detector degrades, while the decorrelating detector [ ] elim- 
inates the user interference by projecting the received signal 
onto the subspace of the signature waveform of each user. The 
decorrelating detector optimizes the near-far resistance among 
linear detectors [ ], although it also amplifies noise. Both 
linear and nonlinear MUDs have sufficiently high complexity 
to preclude their wide adoption in deployed systems. One 
reason is that they both require the number of correlators at 
the receiver front-end to be equal to the number of users in 
the system. 

In an earlier work [ ], we introduced the structure of a 
low complexity reduced-dimension multiuser detection (RD- 
MUD). The RD-MUD exploits the fact that the number of 
active users K is typically much smaller than the total number 
of users N at any given time. Our RD-MUD has a front-end 
that correlates the received signal with M correlating signals, 
with M much smaller than N. The correlating signals are 
formed as linear combinations of the signature waveforms via 
a (possibly complex) coefficient matrix A. Our choice of A 
will be shown to be crucial for performance. The output of 
the RD-MUD front-end can thus be viewed as a projection 
of the MF-bank output onto a lower dimensional detection 
subspace. 

After first developing structures for general linear and 
nonlinear RD-MUDs, we will develop performance analysis 
bounds for two of these structures: the reduced-dimension 
decorrelating (RDD) detector, a linear detector that combines 
subspace projection and thresholding to determine active users 
with a sign detector for data recovery [ ], and the reduced- 
dimension decision-feedback (RDDF), a nonlinear detector 
that combines decision-feedback orthogonal matching pursuit 
(DF-OMP) [4] for active user detection with the sign detector 
for data recovery in an iterative manner. We present theoretical 



probability-of-error performance guarantees for these two 
detectors in terms of the coherence of the matrix A, in a non- 
asymptotic regime with a fixed number of users and active 
users. Our RD-MUD detectors consists of two stages: active 
user detection and data detection of active users. The first 
stage is closely related to [ ]. However, our problem differs 
in that the probability-of-error must consider errors in both 
stages. We derive conditions under which the probability- 
of-error is dominated by errors in the first stage. We do 
not consider optimizing signature waveforms and hence our 
results will be parameterized by the crosscorrelation properties 
of the given signature waveforms. 

The rest of the paper is organized as follows. Section II 
and Section III present the model and the RD-MUD front- 
end, respectively. Section IV introduces the digital detectors 
we propose for RD-MUD. Section V contains the theoretical 
performance guarantees. Section VI contains numerical exam- 
ples, and finally Section VII concludes the paper. 

II. System Model 

Consider a multiuser system with N users. Each user 
is assigned a unique signature waveform from a set 5 = 
{sn{-) : [0, T] ^ R, 1 < n < N}, which are assumed given 
and known, and posess certain properties discussed in more 
detail below. Each user modulates its signature waveform to 
transmit its symbols. The symbols carry information. The 
duration of the signature waveforms T is referred to as 
the symbol time. Define the inner product (or crosscorre- 
lation) between two real analog signals x{t) and y{t) as 
{x{t)^y{t)) = x{t)y{t)dt. The crosscorrelations of 

the signature waveforms are characterized by the Gram matrix 
G, defined as 

[G]ni = {sn{t),si{t)), 1 < n < TV, 1 < / < TV. (1) 
For convenience, we assume that Sn{t) has unit energy: 

IknWiP = {Sn{t),Sn{t)) = 1 for all Tl SO that [G]nn = L 

We also assume that the signature waveforms are linearly 
independent. Hence G is invertible. We consider the syn- 
chronous MUD model that uses Binary Phase Shift Keying 
(BPSK) modulation [ ]. There are K active users with index 
set n e X. The complement set contains indices of all non- 
active users. The symbol of the user n is bn G {1,-1}, for 
n e X. Define a gain factor rn for each user which captures 
the transmitting power and channel gain. We assume is real 
and known to the receiver. The nonactive users can be viewed 
as transmitting with zero power, or equivalently transmitting 
zeros: bn = 0, for n e X^. The received signal y{t) is a 
superposition of the transmitted signals from the active users, 
plus white Gaussian noise w{t) with zero-mean and variance 
a^: 

N 

y{t) = ^rnbnSn{t)^w{t), t G [0,T], (2) 

n=l 

with bn G {1, —1}, n e X, and 6^ = 0, n G X^. The goal of 
multiuser detection (MUD) is to detect the set of active users 
X and their transmitted symbols {bn ' n eX}. In practice the 
number of active users K is typically much smaller than the 
total number of users N, which is a form of user sparsity. 



Therefore, the received signal y{t) consists of only a few 
transmissions from active users. As we show, this user sparsity 
enables us to reduce the number of correlators at the front- 
end and still be able to achieve performance similar to that of 
a conventional MUD using a bank of MFs. To simplify the 
detection algorithm, we assume that K is known. The problem 
of estimating K can be treated separately [ ]. 



y(t) 



Fig. 1: Front-end of RD-MUD. 

III. RD-MUD Front-End 

The RD-MUD front-end, illustrated in Fig. 1, correlates 
the received signal y{t) with a set of correlating signals 
hm{t), m = 1,---M, where M is typically much smaller 
than N. This is in contrast to the conventional matched 
filter (MF) bank, which correlates the received signal with 
the full set of N signature waveforms [ ]. The front-end 
output is processed by either a linear or nonlinear detector 
to detect active users and their symbols, as shown in Fig. 
2 for both linear and nonlinear detectors. The design of the 
correlating signals hm{t) is the key for RD-MUD to reduce 
the number of correlators. To construct these signals, we rely 
on biorthogonal waveforms [ ]. The biorthogonal signals with 
respect to {sn{t)} are defined as a linear combination of 
all signature waveforms using a weighting coefficient matrix 
G-': Sn{t) = Eti[^~\Mt)^ 1 < n < TV, where 
[X]nm denotes the element of a matrix X at its nth row 
and the mth column. Also denote by [x]^ the nth entry of 
a vector x. The biorthogonal signals have the property that 
{sn{t)^ Sni{t)) = Sn,m, for all n, m. The delta function Sn,m 
is equal to one when n = m, and is equal to zero otherwise. 
The correlating signals hni{t) are linear combinations of the 
biorthogonal waveforms with (possibly complex) weighting 
coefficients amn that we choose: 
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hm{t) ^ amnSn{t), 1 < m < M. 



(3) 
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with [Ajmn = dmn 
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- K^ln, • • • , CLMnl ' 



Define a coefficient matrix A G 
and denote the nth column of A as sifi 
n = I,-- - ,A^. The notation denotes the transpose of 
a vector or matrix. We normalize the columns of A so that 
ll^nlP = Em=i Km^nm = 1, whcrc X* is the Conjugate of a 
scalar x. The design of the correlating signals is equivalent to 
the design of the coefficient matrix A for a given {sn{t)}. We 
will use coherence as a measure of the quality of A, which is 
defined as [ ]: 

(4) 



/i = max la^a^ I 



As we will show later, it is desirable that the columns of A 
have small correlation such that fi is small. The output of 



the mth correlator is given by ym = {hm{t)^y{t)). Denoting 
y = [yir ' ' J Vm]'^ , we can derive the output of the RD-MUD 
front-end as (detailed derivations can be found in [2]): 

y = ARb + w, (5) 

where w is a Gaussian random vector with zero mean and 
covariance cr^AG~^A^, R is a diagonal matrix with r^n 
on the diagonal, and b = ,6Ar]^. The notation 

denotes the conjugate transpose of a matrix X. The vector y 
can be viewed as a linear projection of the MF-bank front-end 
output onto a lower dimensional subspace which we call the 
detection subspace. Since there are at most K active users, b 
has at most K non-zero entries. The idea of RD-MUD is that 
when the original signal vector b is sparse, with proper choice 
of the matrix A, the detection performance for b based on y of 
(5) in the detection subspace can be similar to the performance 
based on the output of the MF-bank front-end. 
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Fig. 2: The diagram of (upper) linear detector, and (lower) 
nonlinear detector. 

IV. RD-MUD Detectors 

We now discuss how to recover b from the RD-MUD front- 
end output y of (5) using digital detectors. The model for the 
output (5) of the RD-MUD front-end has a similar form to 
the observation model in the compressed sensing literature 
[7] [5], except that the noise in the RD-MUD front-end output 
is colored due to match filtering at the front-end. Hence, to 
recover b, we can adopt the ideas developed in the context 
of compressed sensing, and combine them with techniques of 
MF-bank detection. 

The linear detector for RD-MUD first recovers active users 
X using support recovery techniques from compressed sensing 
(e.g., [ ]). Given an index set X, Xx denotes the submatrix 
formed by the columns of a matrix X indexed by X, and xj 
denotes the subvector formed by the entries indexed by X. 
Based on the recovered index set of active users X, we can 
write the RD-MUD front-end output model (5) as 



w. 



(6) 



Once the active users are detected, their symbols b^ can be 
detected from (6). This is done by applying a linear transform 
to the front-end output and detecting symbols separately. The 
nonlinear detector for RD-MUD detects active users and their 
symbols jointly (and/or iteratively). 

We will focus on recovery based on two algorithms: (1) 
the reduced-dimension decorrelating (RDD) detector, a linear 
detector that uses subspace projection along with thresh- 
olding [ ] to determine active users and sign detection for 



data recovery; (2) the reduced-dimension decision feedback 
(RDDF) detector, a nonlinear detector that combines decision- 
feedback orthogonal matching pursuit (DF-OMP) for active 
user detection and sign detection for data recovery. DF-OMP 
differs from the conventional OMP [ ] in that in each iteration, 
the binary- valued detected symbols, rather than the real- valued 
estimates, are subtracted from the received signal to form the 
residual used by the next iteration. The residual consists of the 
remaining undetected active users. By subtracting interference 
from the strongest active user we make it easier to detect the 
remaining active users. 

A. Reduced-dimension decorrelating (RDD) detector 

The RDD detector works as follows. As per (5), the front- 
end of the RD-MUD projects the received signal y{t) onto 
the detection subspace as a vector y. By considering the RD- 
MUD front-end output when the input signal is Sn{t), we can 
show that the column of A corresponds to the nth signature 
waveform vector in the detection subspace. Considering the 
detection method of the conventional MUD, a natural strategy 
for RD-MUD is to match the received signal vector y and the 
nth signature waveform vector in the detection subspace by 
computing their inner product, which is given by a^y, n = 
1, • • • , A^. To detect active users, we can rank the magnitudes 
of these inner products and detect the index of the K largest 
as active users: 

I={n: if|K[a^y]| 

is among the K largest of |3?[a^y] |, n = 1, • • • , TV}, 

where denotes the real part of a number x. To detect 
their symbols, we use sign detection: 



(7) 



sgn (rn5R[a^y]) 
0, 



n G X; 
n ^ X. 



(8) 



where sgn(x) denotes the sign of a number x. In detecting 
active users (7) and their symbols (8), we take real parts of 
the inner products because the imaginary part of a^y contains 
only noise and interference. The complexity-per-bit for data 
detection of the RDD detector is proportional to M. Since 
M < TV in RD-MUD, the complexity for data detection of the 
RDD detector is on the same order as that of the conventional 
linear MUD detector. But the RDD detector requires much 
lower decorrelation complexity than the conventional linear 
detector. 

B. Reduced-dimension decision feedback (RDDF) detector 

The RDDF detector detects active users and symbols it- 
eratively. It starts with an empty set as the initial estimate 
for the set of active user X^, zeros as the estimated symbol 
vector b*^^^ = 0, and the front-end output as the residual vector 
v*^^^ = y. Subsequently, in each iteration /c = 1, i^T, the 
algorithm selects the column an that is most highly correlated 
with the residual v*^^~^^ as the detected active user in the /cth 
iteration, with the active user index: 



n/e = argmax 



.(fe-i)i 



(9) 



This index is then added to the active user set X^^^ = X^^ U 
{n/c}. The symbol for user Uk is detected with other detected 



symbols staying the same: 

sgn(SR[r„ a^v(fc-i 



^(fc-l) 



n Uk] 
n ^ Uk. 

Then the residual vector is updated through 



(10) 



(11) 



The iteration repeats K times (we show in [ ] that with high 
probability DF-OMP never detects the same active user twice), 
and finally the active user set is given by X = X^^^ with the 
symbol vector bn = n = 1, • • • ^N. The complexity- 

per-bit of the RDDF detector is proportional to MK. Since 
M < N, this implies that the complexity for data detection 
of the RDDF detector is on the same order as that of the 
conventional DF detector. But the RDDF detector requires 
much lower decorrelation complexity than the conventional 
DF detector. 

C. Reduced-Dimension MMSE (RD-MMSE) Detector 

Similar to the MMSE detector of the conventional MUD, a 
linear detector based on the MMSE criterion can be derived 
for the reduced-dimension model (6) (see [8] for derivations). 
The RD-MMSE detector detects the set of active users X first 
by a support recovery method and then detects symbols as: 

I sgn([%Af (A^R| Af + a^AG-^A^j-iyln), net- 



0, 



(12) 



D. Maximum likelihood detector 



The optimal detector that minimizes the probability-of- 
error for the RD-MUD output is the nonlinear maximum 
likelihood detector. The maximum likelihood detector finds 
the active users and symbols by minimizing the likelihood 
function, or equivalently, minimizing the quadratic function 
||(AG"^A^)-i/2(y _ ARb)|p. This is also equivalent to 
solving the following integer optimization problem 



max 2y^ (AG" ^ A^ ) " ^ ARb 

6^G{-1,0,1} 

- b^RA^(AG"^A^)-^ARb, 



(13) 



where 6^ = corresponds to the nth user being inactive. 

E. Choice of A 

The coefficient matrix A is our design parameter. In Section 
IV- A and Section IV-B we have shown that both the RDD and 
RDDF detectors are based on the inner products between the 
projected received signal vector and the columns of A. Hence, 
intuitively, for the RDD and RDDF detectors to work well, 
the inner products between columns of A, or its coherence 
defined in (4) should be small. In the following we consider 
the random partial discrete Fourier transform (DFT) matrix, 
whose coherence is small and it is formed by randomly 
selecting rows of a DFT matrix F: [Fj^m = 
normalizing the columns of the sub-matrix, where i 



nm 



V. Performance of RD-MUD 

In the following, we present conditions under which the 
RDD and RDDF detectors can successfully recover active 



users and their symbols. The conditions depend on A through 
its coherence and are parameterized by the crosscorrelations 
of the signature waveform through the properties of the matrix 
G. Our performance measure is the probability-of-error, which 
is defined as the chance of the event that the set of active users 
is detected incorrectly, or any of their symbols are detected 
incorrectly: 

Pe = PiX + r)-r P{{X = X} n {b ^ b}). (14) 

We will show that the second term of (14) is dominated by the 
first term. The noise plays two roles in the Pg of (14). First, the 
noise can be sufficiently large relative to the weakest signal 
such that a nonactive user is determined as active; second, 
the noise can be sufficiently large such that the transmitted 
symbol plus noise is detected in an incorrect decision region 
and hence decoded in error. The first error term in (14) is 
related to the probability-of-error for support recovery (see, 
e.g. [9]). There are two major differences in our results on 
this aspect of RD-MUD performance relative to those previous 
works. First, although noise in the analog signal model (2) is 
white, matched filtering at the RD-MUD front-end introduces 
colored noise in (5). Second, we take into account the second 
term in (14), which has not been considered in previous work. 
Define the largest and smallest channel gains as 



I I A I I 

Tmaxl max|rn|, 

n—1 



A ^ 

= min \rn 

n—l 



(15) 



Our main result is the following theorem: 



Theorem 1. Let b G 



pATxl 



be an unknown deterministic 



symbol, bn G { — 1,1}, n G X, and bn = 0, n G X^, n = 
1, • • • , N. Assume that the number of active users K is known. 
Given the RD-MUD front-end output y = ARb + w, where 
A G ^ ^ and G G ^ ^ are known, and w is a Gaussian 
random vector with zero mean and covariance cr'^AG~^A^ , if 
the columns of A are linearly independent and the coherence 
of A (4) satisfies the following condition: 



{2K - 1)m| ' max 



(16) 



for some constant a > 0, and N *^^+"^[7r(l 
a)log7V]~^/^ < 1, where 



T ^ (7^2(1 + a) log TV • ^ A„,ax(G-^) • ^max (a^AA^a^), 

(17) 

then the probability-of-error (14) for the RDD detector is 
upper bounded as: 

Pe < 7V""[7r(l + a)log7V]-^/^ (18) 

If the columns of A are linearly independent and the coherence 
of A (4) satisfies a weaker condition: 



(2i^-l)/i|r^in| >2r, 



(19) 



for some constant a > 0, and A/'~^^+"^ [7r(l + 
a)logA^]~^/^ < 1, then the probability-of-error (14) for the 
RDDF detector is upper bounded by the right hand side of 
(18). 

The proof for Theorem 1 is given in [8]. The key idea of the 
proof is to find a uniform bound for the tail probability of the 



correlator output noise. Note in Theorem 1 that the condition 
of having small probability-of-error for the RDDF detector is 
weaker than for the RDD detector. Based on the coherence of 
the random partial DFT matrix, we can prove the following 
corollary to Theorem 1 (see [ ] for more details): 

Corollary 1. Consider the setting of Theorem 1, where A 
is a random partial DFT matrix A. Suppose the number of 
correlators satisfies the following lower bound for the RDD 
detector 

2 



N = 100, K = 2 



M > 4 



(2i^-l)|r„ 



2t 



(21ogiV + c), 



(20) 



or satisfies the following smaller lower bound for the RDDF 
detector 

2 



M > 4 



(2i^-l)|rn 



2r 



(2 log TV + c), (21) 



for some constants c > and a > 0, and |rmin| > 2r for r 
defined in (17), then the probability-of-error Pg of the RDD 
detector or the RDDF detector is bounded by 

1 - (1 - 7V-"[7r(l + a) log A/']-^/^)(l - 26"^), (22) 

for some constant a > 0. 

This corollary says that to attain a small probability-of- 
error, the number of correlators needed by the RDD and 
RDDF detectors is on the order of logA^, which is much 
smaller than that required by the conventional MUD using 
the MF-bank, which is on the order of TV. 

VI. Numerical Examples 

As an illustration of the performance of RD-MUD, we 
present an numerical example using the RDD detector. The 
results are obtained from 10^ Monte Carlo trials. For each 
trial, we generate a Gaussian random noise vector w as well 
as a random partial DFT matrix for A, and form the signal 
vector according to (5). To simplify, we assume that the gains 
for all the users are the same! |rniin| — |^max| — 1- Assume 
the signature waveforms are orthogonal (G = I). In this case, 
the noise in (5) is white. We define the signal-to-noise-ratio 
(SNR) as |rmin|V^^ = 1/^^- We also assume TV = 100 and 
K = 2. Fig. 3 shows Pg versus M for the RDD detector as 
SNR increases. The counterpart of RD-MUD with the RDD 
detector in the conventional MUD setting is the decorrelating 
detector (when no subspace projection happens, i.e., if we 
let A = I in (5)). For each SNR, as M increases, the Pg 
of the RDD detector approximates that of the conventional 
decorrelating detector. Also with higher SNR, the Pg of the 
RDD detector decreases faster with increasing M. When SNR 
is sufficiently high, the number of correlators required by the 
RDD detector to achieve a small Pg is much fewer than N. 
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Fig. 3: Performance of the RDD detector, Pg versus M for 
different SNRs, when the signature waveforms are orthogonal, 
i.e., G = I. The dashed lines show Pg for the conventional 
decorrelating detectors at the corresponding SNR. When SNR 
is greater than 15 dB, (with = 100 correlators) the 
probability-of-error of the decorrelating detector is less than 
10-^ 



method reduces the number of correlators at the front-end of a 
MUD receiver by exploiting the fact that the number of active 
users is typically much smaller than the total number of users 
in the system. Motivated by the idea of analog compressed 
sensing, the RD-MUD front-end projects the received signal 
onto a lower dimensional detection subspace by correlating 
with a set of correlating signals. We proved that when the 
random partial DFT matrix is used to construct correlating 
signals for RD-MUD, the number of correlators is on the 
order of log of the number of users in the system, which 
is much smaller than that required by the conventional MUD. 
Numerical examples validated our theoretical results. 
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VII. Conclusions 

We have developed families of digital detectors for 
the reduced-dimension multiuser detection (RD-MUD), and 
proved performance guarantees for two specific detectors: 
the reduced-dimension decorrelating (RDD) detector and the 
reduced-dimension decision feedback (RDDF) detector. This 



